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Multi-particle entanglement is a fundamental 
feature of quantum mechanics |lj and an essential 
resource for quantum information processing |2j 
and interferometry [3, 4|. Yet, our understanding 
of its structure is still in its infancy. A system- 
atic classification of multi-particle entanglement 
is provided by the study of equivalence of en- 
tangled states under stochastic local operations 
and classical communication [5J [6]. Determining 
the precise entanglement class of a state in the 
laboratory, however, is impractical as it requires 
measuring a number of parameters exponential in 
the particle number. Here, we present a solution 
to the challenge of classifying multi-particle en- 
tanglement in a way that is both experimentally 
feasible and systematic, i.e., applicable to arbi- 
trary quantum systems. This is achieved by asso- 
ciating to each class an entanglement polytope — 
the collection of eigenvalues of the one-body re- 
duced density matrices of all states contained in 
the class. Determining whether the eigenvalues 
of an entangled state belong to a given entangle- 
ment polytope provides a new criterion for multi- 
particle entanglement. It is decidable from a lin- 
ear number of locally accessible parameters and 
is robust to experimental noise. We describe an 
algorithm for computing entanglement polytopes 
for any number of particles, both distinguishable 
and indistinguishable. Furthermore, we illustrate 
the power of entanglement polytopes for witness- 
ing genuine multipartite entanglement and relate 
them to entanglement distillation. The polytopes 
for experimentally relevant systems comprised of 
either several qubits or bosonic two-level systems 
are explained. 

It is an oft-repeated piece of folklore that when it comes 
to the relation between a global quantum state and its 
local reductions, "the whole is greater than the sum of 
its parts" [7, ,8] . This belief seems particularly relevant 
to entangled many-body quantum systems, routinely con- 
trolled experimentally [3J l9Tlll|. where the number of 
global parameters is exponentially larger than number 
of locally accessible ones. We challenge this intuition by 
showing that single-particle information alone can serve 
as a powerful witness to multi-particle entanglement. Our 
methods apply to states that are sufficiently pure — a fact 
that itself can be efficiently verified experimentally. 

Formally, a pure state is said to be entangled if it cannot 
be written as a product \ip) ^ \ip^) ® . . . ® |V> (Ar) ) 0- 



As entanglement theory explores the non-classical aspects 
of quantum states, it is natural to think of two states as 
belonging to the same entanglement class if they can be 
converted into each other with finite probability of success 
using only local operations and classical communication 
(stochastic LOCC, or SLOCC) 0(5]. 

For small systems, these entanglement classes are well- 
understood. In the simplest scenario of three qubits 
(two-level systems) , there exist two classes of genuinely 
entangled states of very different nature: the first contains 
the famous Greenberger-Horne-Zeilinger (GHZ) state 
^g(lttt) + 144-4-))) which exhibits a particularly strong 
form of quantum correlations [lj; the second contains 
the W state ^(Ittl) + It4t) + I4tt)) E- Whereas a 
single GHZ state can be approximately converted into 
a W state, the converse is not true — implying stronger 
entanglement of the former. Already for four subsystems 
there exist infinitely many entanglement classes |12| . and 
the number of parameters required to determine the class 
grows exponentially with the number of particles. As a 
result, only sporadic facts are known for larger systems, 
despite the enormous amount of literature dedicated to 
the problem. 

The interest in entanglement classes stems not only 
from their role in quantum information. Indeed, connec- 
tions to relativistic transformations and even to black 
holes have been found [13]. Mathematically, the problem 
can be formulated in terms of invariant theory, studied 
since the 19th century — Cayley's hyperdeterminant, e.g., 
appears as the 3-tangle [14J . Similar techniques underpin 
modern developments, such as the geometric approach to 
the P vs. NP problem |15j . 

Given the natural interpretation and the broad appli- 
cability of the concept of SLOCC equivalence, it is highly 
relevant to develop a more tractable approach to its study. 
The challenge is to identify a method that (i) involves only 
a manageable number of parameters, (ii) is experimen- 
tally feasible, (in) can systematically be used for systems 
of arbitrary size, and (iv) produces succinct results. 

In this work, we propose a solution that meets all 
the above desiderata. It is achieved by establishing a 
connection between the entanglement problem and a 
seemingly unrelated question: the one-body quantum 
marginal problem (or N -representability problem in chem- 
istry [TB]). This fundamental problem about quantum cor- 
relations asks which set of single-particle density matrices 
(p^\ . . . , p^ N ') can appear as the set of reduced density 
matrices of a globally pure quantum state. Its solution is 
easily seen to depend only on the eigenvalues of the 
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Figure 1. Entanglement polytopes as witnesses, (a) An 

entanglement polytope contains all possible local eigenvalues of 
states in the entanglement class (the W class and its polytope 
are shown in dark gray), (b) For a sufficiently pure quantum 
state p prepared in the laboratory, local tomography is per- 
formed to determine its local eigenvalues, (c) The indicated 
eigenvalues are not compatible with the W class, hence p must 
be entangled of GHZ type. 

densities and allows for an extremely elegant mathemati- 
cal description: the set of possible (A^, . . . , A^) forms 
a convex polytope |17H19| whose defining inequalities can 
be computed algorithmically |181 ITU] . (For fermions, the 
most famous such inequality is nothing but the Pauli 
principle |20|). 

The crucial observation is that the set of allowed eigen- 
values may depend on the entanglement class of the global 
state. Thus, we propose to associate to every class the col- 
lection of eigenvalues of the one-body reduced density ma- 
trices of the states in its closure. We show that these sets 
also form convex polytopes, which we call entanglement 
polytopes. Entanglement polytopes immediately lead to a 
local criterion for witnessing global multi-particle entan- 
glement: If the collection of eigenvalues (A^, . . . , A^) of 
the one-body reduced density matrices of a pure quantum 
state does not lie in an entanglement polytope Ac, 
then the given state cannot belong to the corresponding 
entanglement class C (Figure 1): 

(A«...,AW)£ A c => \^)iC. (1) 

Phrased differently, the criterion allows us to witness 
the presence of a highly entangled state by showing that 
its eigenvalues are incompatible with all less-entangled 
classes. The conclusions drawn from (|T|) are the strongest 
possible statements about the entanglement class of the 
state based on local information alone. 

Crucially, there are always only finitely many entangle- 
ment polytopes and we describe a method to compute 
them. They naturally form a hierarchy: if a state in 
the class C can be converted into a state from T> then 
Ac 3 Ax> . This reflects the fact that states of the first 
class are more powerful for quantum information process- 
ing. We illustrate the method by considering relevant 
examples (cf. Supplementary Information): 

3 qubits |21| : There are two full-dimensional entan- 
glement polytopes: one corresponding to the W class 
(the dark upper pyramid in Figure 1) and one to the 




Figure 2. Cuts through two entanglement polytopes for 
four qubits. Each row represents a single four-dimensional 
entanglement polytope. The coordinate A^lx is fixed in every 
column — to values 0.5, 0.66, 0.83, and 1, respectively — and the 
remaining three variables are shown. The first row corresponds 
to the class L a4 for a = in the characterization of |12| . The 
second row is the four-partite W state — one can clearly identify 
the characteristic "upper-pyramid form" as explained in the 
text. Several properties that previously had to be computed 
algebraically (e.g., in |12j ) can be read off directly: E.g., the 
final column corresponds to the situation when the first system 
has been projected onto a generic product state. In this case, 
we recover a state vector on the remaining three sites, which 
can be seen from its polytope to be of GHZ type in the first 
row, and of W type in the second row. 

GHZ class (the entire polytope, i.e., the union of both 
pyramids). The tip of the upper pyramid constitutes 
a polytope by itself, indicating a product state. Three 
further one-dimensional polytopes are given by the edges 
emanating from this vertex. They correspond to the three 
possibilities of embedding a Bell state and a single-qubit 
state into three systems. 

Thus, eigenvalues in the interior of the polytope are 
compatible only with W and GHZ classes, i.e., genuine 
three-partite entanglement. Likewise, if the eigenvalues 
are contained in the lower pyramid 

\nax "F -^rnax + ^LL < 2, (2) 

then by ([TJ the state cannot be contained in the W class — 
we have witnessed GHZ-type entanglement. 

4 qubits: There exist 9 infinite families of entanglement 
classes, each described by up to four complex parameters 
|12j that are not directly accessible. Hence, arguably, 
a complete classification is too detailed to be practical. 
In contrast, the polytope method strikes an attractive 
balance between coarse-graining and preserving structure 
[22]: 

Up to permutations, there are 13 entanglement poly- 
topes, 7 of which are full-dimensional and correspond 
to distinct types of genuine four-partite entanglement. 
One example is given by the 4-qubit W class, which, 
in complete analogy to the previous case, consists of 
an "upper pyramid" of those eigenvalues that fulfill 

A^L + a£L + ASL + ASL > 3 (Figure 2). 

N qubits & genuine multipartite entanglement: \tp) is 
bi-separable if it factorizes with respect to some partition 
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Figure 3. Witnessing genuine TV-partite entanglement. 

The figure displays the two-dimensional cut through the six- 
qubit eigenvalue polytope where we fix A^ n = . . . = A^ ; ' n = 
1/8. If the collection of local eigenvalues does not belong to the 
union of the biseparable entanglement polytopes (dark gray 
region) then any corresponding pure quantum state necessarily 
contains genuine six-qubit entanglement. 



of the qubits into two sets. Otherwise, \tp) is genuinely N- 
partite entangled. Entanglement polytopes of bi-separable 
states do not account for all possible collections of eigen- 
values. Hence, verifying that the local eigenvalues of a 
state do not lie in any bi-separable polytope provides a 
criterion for genuine TV-partite entanglement (Figure 3). 

N bosonic qubits: In the case of two-level systems in a 
Bose-Einstein condensate (e.g. [33]), all reduced density 
matrices coincide and are described by a single number: 
the maximal local eigenvalue. The entanglement poly- 
topes are intervals from suitable constants jc to the value 
one. The most entangled (see below) states in any class 
are given by the Dicke basis states (7c > 0.5) and they 
include the GHZ state for 7c = 0.5. 

Entanglement polytopes also give rise to single-number 
entanglement measures inferred from local information. 
Witnessing genuine multipartite entanglement as indi- 
cated above serves as one example. For a second exam- 
ple, consider the linear entropy of entanglement E(tp) = 

1 - ^Ei/li tr (j) ) 2 used > e -S-> in metrology [H]. The 
polytope method allows us to bound the maximal entropy 
of entanglement distillable from the global state, even 
though only an exponentially small fraction of the state 
is known. For this, we note that E(tp) corresponds to 
the Euclidean length of the vector of local eigenvalues: 
shorter vectors imply more entanglement. In particular, 
the quantum states of maximal entropy of entanglement 
in a class C map to the point of minimal distance to the 
origin in the entanglement polytope Ac- Therefore: If 
the local eigenvalues of a given state lie only in entan- 
glement polytopes that come close to the origin, then a 
high amount of entanglement can be distilled from the 
state. This observation can easily be turned into a more 
quantitative solution (see Supplementary Information). 

Generalizing an algorithm described in [25] , we even ob- 
tain a physical scheme for a distillation procedure without 
the need to learn the global state: Repeatedly perform 
local tomography and apply the stochastic local operation 




Figure 4. Entanglement distillation for a quantum 
state of GHZ class. By iterating the distillation procedure, 
a quantum state with maximal linear entropy of entanglement 
is reached asymptotically. 

corresponding to the local density matrices themselves. 
If the operation has been performed successfully, the en- 
tropy of entanglement increases in every step and the 
eigenvalue vector flows to the closest point to origin in its 
polytope (Figure 4). 

The described methods are efficient in the number 
of parameters to be measured, as well as robust to ex- 
perimental noise and can readily be realized in current 
experiments. Since quantum states p in the laboratory 
are always mixed, one first estimates the purity trp 2 , 
e.g., by using a "swap test", which only requires two-body 
measurements [26 . In a second step, one performs local 
tomography of the one-body density matrices. The total 
number of measurements is proportional to the number 
of particles N. If the purity is close to unity and the 
local eigenvalues sufficiently far from the boundary of any 
entanglement polytope, then two rigorous statements can 
be made: First, the density matrix p has high fidelity 
with a pure state \ip) that has the entanglement properties 
determined by the polytope method. Second, the mixed 
state p cannot be written as a convex combination of pure 
states in classes excluded by the polytope analysis. In 
these ways, all statements can be lifted to the case with 
experimental noise. 



METHODS SUMMARY 

We now list the methods used to prove the claims above. 
Fully worked out examples, proofs of all technical claims, 
as well as the mentioned algorithms can be found in the 
Supplementary Information. 

We consider an entanglement class C and its correspond- 
ing entanglement polytope Ac . In order to see that Ac is 
indeed a convex polytope, we employ the description of 
SLOCC operations as elements A\ ® . . ,®Ajv of the group 
G of invertible local operators [BJ. By using results from 
algebraic geometry, one can express Ac in terms of the 
decomposition into irreducible G-representations of the 
ring of polynomial functions on the closure C, or, equiva- 
lently, in terms of covariants |27] . Intriguingly, this can be 
understood as the semi-classical limit of a geometric quan- 
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tization in the sense of Kostant and Souriau. Covariants, 
which generalize the notion of an invariant [28_ , are vector- 
valued polynomial functions on the Hilbert space % which 
transform like an irreducible representation of the group 
G. The statement that Ac is a convex polytope follows 
from the fact that covariants form a finitely-generated al- 
gebra [57] • We show that, crucially, finite generation also 
implies that the entanglement polytopes can be computed 
algorithmically by using techniques from computational 
invariant theory [29J . For the purpose of computing exam- 
ples, we use known lists of covariants and the description 
of the full marginal polytope for qubit systems. 



For the entanglement distillation procedure, we show 
that E(p) can be increased by following its gradient flow, 
which is given by the infinitesimal action of the collec- 
tion of reduced density matrices. Any local maximum 
attained in this way is automatically a global maximum. 
Mathematically, the limit points are critical points in the 
canonical Kempf-Hesselink-Kirwan stratum with respect 
to the equivariant Morse gradient flow |30j . They cor- 
respond geometrically to the point in the entanglement 
polytope which has minimal Euclidean distance to the 
origin. In practice, the flow is approximated by finite time 
steps, each of which corresponds to an SLOCC operation. 
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SUPPLEMENTARY INFORMATION 



Introduction 



In the following, we will give rigorous proofs of the main properties of entanglement polytopes ( § 3 ) and describe 



an algorithmic method for their computation, which we illustrate with several worked examples (§5 1. Moreover, 



we elaborate on the properties of the linear entropy of entanglement and derive the gradient flow procedure for 



entanglement distillation that has been sketched in the main body of the article (§4|. Our main technical tools are 
Brion's invariant-theoretic description of moment polytopes [T] seen through the lens of non-reductive group actions 
[2 , and Kirwan's analysis of the equivariant Morse gradient flow for the norm square of a moment map [3J. 

For clarity of exposition, we will describe our results for a system composed of N distinguishable subsystems with 
dk degrees of freedom, respectively, and Hilbert space 

% = C dl ® . .. <g> C dN . 

In the context of multi-particle entanglement, we will think of each of the N subsystems as corresponding to an 
individual particle. However, the subsystems can be of more general nature and, e.g., describe different degrees of 
freedom such as position and spin. All results can be adapted to systems composed of bosons or fermions by replacing 



% with the (anti) symmetric subspace (§5 1. We shall denote by 

p(H) = {/> = |Wh WW = i} 

the projective space of pure states; the expectation value of an observable O in a pure state p is given by tr(pO) = 
(tp\0\ip)- The effective state of the fc-th particle is described by the (one-body) reduced density matrix p^ k \ which by 
definition reproduces the expectation values of all one-body observables 0^ k \ 

tv(p^O^) = tr(p (l®^" 1 ) (g) (k) <g> = (^l®**- 1 ) ® {k) ® l®( JV - fc )|^) (A.l) 

The one-body reduced density matrices p^ k ' are in general mixed states, i.e., convex mixtures of pure states. Formally, 
they are positive semidefinite Hermitian operators of trace one. By the spectral theorem, each can be diagonalized, 
and the local eigenvalues = X(p^) are vectors of non- negative numbers summing to one (which we assume by 
convention to be weakly decreasing). 



2. Multi-Particle Entanglement and Stochastic Local Operations and Classical Communication (SLOCC) 

A pure quantum state p — of a multi-particle system is called entangled if it cannot be written as a tensor 

product [3] 

It is easy to see that p is unentangled, or separable, if and only if all its one-body reduced density matrices p( fe ' are 
pure states, that is, if and only if A^ fe ) = (1, 0, . . . , 0) for k = 1, . . . , N. 

In order to classify the entanglement present in a given multipartite quantum state p, it is useful to compare its 
capability for quantum information processing tasks with that of other quantum states. Specifically, we shall think of p 
to be at least as entangled as any other quantum state p' that can be produced from a single copy of p by performing 
a sequence of stochastic local operations (i.e., local operations which succeed with some positive probability, e.g., 
by post-selecting on a certain measurement outcome) and classical communication (SLOCC) [5J|B]. If, conversely, 
p can also be produced from p' then we can think of the two states as possessing the same kind of multi-particle 
entanglement. In this way the set of quantum states is partitioned into equivalence classes. For pure states p = IV'KV'I 
and p' = \ip')(ip'\, it has been shown that they are equivalent under SLOCC if and only if there exist invertible 
operators such that 

=(A«<8...®AW)M. 

Indeed, the operators A^ 5 ) can be defined by following a successful branch of a conversion protocol. Conversely, given 
operators A^ it suffices to perform successive local POVM measurements with Kraus operators 

S (k) = y/tfK)A< k \ F« = ^l-p( fc )(AW)U( & ), (A.2) 
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where p^ is a suitable normalization constant. If all measurements succeed (outcomes S^) then p' 

has been successfully distilled from p [B]. Mathematically, it is convenient to define an action of the Lie group 
G = SL(d 1 ) x ... x SL(d/v) of invertible local operators of determinant one on the projective space ¥(W) by 

5 ■ lwl := w^^mw ' ( } 

where g = (g^\ . . . , g^) G G. Then the result of [6] states that two pure states p and p' are equivalent under SLOCC 
if and only if G ■ p = G ■ p' . In other words, the SLOCC entanglement class containing p is just the orbit C — G ■ p. 
Clearly, the unentangled states form a single SLOCC class. 

The closure C = G ■ p of an entanglement class contains in addition those quantum states which can be arbitrarily 
well approximated by a state in the class. In this way, the closure of an entanglement class can be given a similar 
operational interpretation as the class itself. While the entanglement classes partition the set of multi-particle quantum 
states, their closures naturally form a hierarchy. This is because they are stable under SLOCC operations; indeed, it is 
immediate that p' G G ■ p implies G • p' C G • p. In particular, every entanglement class contains in its closure the 
class of unentangled states. 

In summary, stochastic local operations and classical communication provide a systematic way of studying multi- 
particle entanglement. However, it is immediate from the fact that the dimension of G grows only linearly with N that 
there is generically an infinite number of distinct SLOCC entanglement classes labeled by an exponential number of 



continuous parameters (cf . § 5 ) . It is therefore necessary to coarsen the classification in a systematic way in order to 
arrive at a tractable way of witnessing multi-particle entanglement. 

We note that an extraordinary amount of research has been devoted to the task of classifying SLOCC entanglement 
and identifying it experimentally. The field is far too large to allow for an exhaustive bibliography. For reviews on the 
general theory, see [TJ [5]; a review focussing on detection is Methods from classical invariant theory have long been 
used to analyze entanglement classes, see, e.g., |10fE5] and references therein. 



3. Entanglement Polytopes 

a. Definition 

The entanglement polytope of an entanglement class C is by definition 

A c = {(A«,...,AW) : XW=A(p«) )P GC}, 

the set of tuples of local eigenvalues of the quantum states in the closure of the entanglement class. We will show 
below that this set is in fact a convex polytope. The set of entanglement polytopes forms a hierarchy which coarsens 
the hierarchy of the closures of entanglement classes described above. Namely, 

C C V => A c C A -p. 

In particular, if a pure quantum state p = \ip){ij)\ is contained in an entanglement class C then its collection of local 
eigenvalues is contained in the corresponding entanglement polytope Ac- That is, 

peC (A« ...,AW)e A c . 

Equivalently, if the collection of local eigenvalues is not contained in the entanglement polytope then the quantum state 
is necessarily in a different entanglement class. This establishes (1), our main criterion for witnessing multi-particle 
entanglement. 



b. Invariant- Theoretic Description 



In order to analyze the properties of entanglement polytopes, it is useful to introduce the map 

$:p = |^<VI^(p (1) ,--- ! P W ) 

which assigns to a quantum state the collection of its one-body reduced density matrices. Given a product of local 
unitaries U = C/ (1) ® . .. ® U^ N \ it follows from ((All that 



$(f/pt/t) = (U^p^(U^)\ . . . , C/WpW(f/W)t). 



(A.4) 
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We can therefore jointly diagonalize the reduced density matrices by applying suitable local unitaries. On the other 



hand, we observe that the action p i— > UpU* is simply the restriction of (A. 3) to the subgroup of local unitaries 
K = SU(<ii) x ... x SU(<i/v) of the group G = SL(<ii) x . . . x SL(djv) introduced in § 2 (the denominator is equal to 
unity since any unitary is norm-preserving), so that the entanglement class is left unchanged. We conclude that for 
every quantum state p £ C there exists a quantum state p' in the same entanglement class whose reduced density 
matrices are diagonal in the computational basis. As we may identify diagonal density matrices with their collection of 
eigenvalues (for definiteness, we shall require the diagonal entries to be arranged in weakly decreasing order), each 
entanglement polytope can be written as an intersection 

A c = $(C) n D h (A.5) 

where we denote by D± the set of tuples of diagonal density matrices with weakly decreasing entries. 

In mathematical terms, the map $ can be identified with the canonical moment map for the action of K on the 



projective space ¥(H) pH [55]. Indeed, by (A.4) it is equivariant with respect to the action of K on tuples of reduced 
density matrices, which can be identified with the coadjoint action of K on the dual of its Lie algebra, and its 
components are Hamiltonian functions with respect to the Fubini-Study symplectic form wfSj 

d X) (iY p ) =itv([Y,p]X) =itr(p[X,Y]) = u FS (iX p , iY p ), (A.6) 

where X and Y are Hermitian operators on % and where Z p — jjj\ t _ Q e tz • p — Jj| t _ e tZ {il)\e tZ / ll e * Z IV 7 ) 111 ^ s 
the tangent vector generated by the infinitesimal action of an arbitrary operator Z; in particular, iX p — i[X, p] and 
X p = {X,p} — 2tr(pX)p. Formally, the above equation should be restricted to local observables X of trace zero, since 
only these correspond to elements in the Lie algebra of the group K. 

It follows that ( |A.5[ ) is just the definition of the moment polytope, or Kirwan polytope, for the subvariety C C P('H), 
with corresponding to a choice of positive Weyl chamber. It is a celebrated result in mathematics that the moment 
polytope of a compact symplectic manifold is a convex polytope [57]. For the singular case of G-stable irreducible 
projective subvarieties of projective spaces, which is a category of spaces that includes the closures of any entanglement 
class, this has been established by Brion pQ. 

Brion's proof of convexity relies on an alternative, invariant-theoretic description of the moment polytopes. In 
order to apply his results to the case of entanglement polytopes, we need some notions from representation theory 
and invariant theory. Recall that the finite-dimensional polynomial irreducible representations of the group GL(d) of 
invertible d x d- matrices are labeled by integer vectors p £ Z> with weakly decreasing entries — their highest weight, 
or Young diagram [28 . The irreducible representations of SL(d) can be obtained by restriction, and they will be 
denoted by Vi. Finally, a covariant of degree n and weights pP~\ ■ ■ ■ , pS N ^ is by definition a G-equivariant map 

n = c dl ® . . . c dN — ► v*fa <g> . . . ® v dN N) 

whose components are homogeneous polynomials of degree n. Note that the right-hand side is a G-representation in 

the obvious way. We shall require each p^ to satisfy Y^j = n \ this can always be satisfied and does not impose 
an additional restriction. Moreover, we observe that even though covariants are defined on the Hilbert space %, their 
non-vanishing can be well-defined on subsets of the projective space P('H) by taking any representative vector, since 
they are homogeneous functions; we will do so without further mention. We can now specialize the main result of [1_ 
to our setting: 

Theorem 1. The entanglement polytope of an entanglement class C — G ■ p is equal to 

— (p^, ■ ■ ■ , jttW) : 3 a covariant of degree n and weights pS k ^ which is non-zero at p\ 
n J 



Proof. By the above discussion and (A.5 1, we have seen that is equal to the moment polytope of the G-orbit closure 
C. In pQ it is shown that this moment polytope can be described as the closure of the set of points (jr- 1 ' , ■ ■ ■ , pf- N ') /n 
for which there exists a non-zero G-equivariant morphism 



N) 



which is homogeneous of degree n (i.e., whose components are homogeneous polynomial functions of degree n). Here, 
we denote by C — {\ip} : |V'}(V'I / (V'lV') £ £•} the affine cone over C. Evidently, every covariant which is non-zero at p 
restricts to a map on C which does not vanish completely; this shows one inclusion. For the converse, we first observe 
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that every non-zero map ip as above is automatically non-zero at p, for it is G-equivariant and the affine cone over 
C = G ■ p is an open subset of C. It remains to show that it can be extended to a G-equivariant map on all of %. This 
can be seen as follows by using some basic algebraic geometry and representation theory [2S] [3D] : 

Since the map (p is G-equivariant and its range an irreducible representation, it is completely equivalent to consider 
instead of tp its component <fi = (ip, v) evaluated at a highest weight vector v with respect to some maximal unipotent 
subgroup U C G. The function is a [/-invariant homogeneous polynomial function of degree n, i.e., a [/-invariant 
element of the degree-n piece of the homogeneous coordinate ring 

C[C]„ = C[H] n /In, 

where C[H] denotes the algebra of polynomial functions on T-L, C[H] n the subspace of homogeneous functions of degree 
n, and I„ = {p £ C["H] n : p(C) = 0} the corresponding part of the vanishing ideal of C. By averaging over the maximal 
compact subgroup K C G, we can find an orthogonal complement J„ to /„ that is G-invariant (this is oftentimes 
called Weyl's unitary trick). Therefore, the natural map 

In Q C[7Z]„ -> C[C]„ 

is a G- isomorphism, and it can be used to find a [/-invariant extension of <fi to all of T-L. By reversing the above 
procedure, we obtain a corresponding covariant which G-equivariantly extends ip to all of %. □ 



c. Properties and Computation 

The set of covariants of fixed degree and weight form a vector space. Moreover, given any two covariants P. Q of 
degree n,m and weights p^ k \ v^ k \ we can form their product by defining 

(P o Q)(|^» - (tt« ® . . . ® ttW) (P(|V)) «8 Q(|V») • (A.7) 

where 7r( fe ) : V^ k) ® Vt^-, — > V^ k)+j y W denotes the projection onto the unique irreducible representation in the tensor 
product of two irreducible representations whose highest weight is the sum of the highest weights of the factors. The 
function P o Q is a covariant of degree n + m and weig hts jl^ + pi k \ It follows that the entanglement polytopes are 
convex, since the point corresponding to P o Q, 

(pP + . . . , m W + = f - ■ ■ ■ , + f 1 ■ ■ • . * N) )) • ( A - 8 ) 

is a convex combination of the points corresponding to the covariants P and Q. If we allow for formal linear combinations 
of covariants of different degree or weight, we thus obtain a graded algebra, and it is well-known that this algebra is 
finitely generated, i.e., there exist finitely many covariants Pi, ... , P m such that any other covariant can be expressed 
as a linear combination of monomials in the Pj. In the literature, this is typically phrased as the well-known statement 



that the algebra C['H]' 7 of [/-invariant polynomials is finitely generated [3T][32]; (A.7) then corresponds to the ordinary 
multiplication of functions. As noted in [T], finite generation implies that moment polytopes are compact convex 
polytopes. Using [Theorem l] we can slightly sharpen this result for the case of entanglement polytopes: 

Corollary 2. The entanglement polytope of an entanglement class C — G ■ p is given by the convex hull 



Ar =< 



I(#\...^):P j( ,)*0 



where we denote by nj the degree and by p^ the weights of a generator Pj (j — 1, . . . , m). In particular, Ac is a 
compact convex polytope. 

Proof. Consider a monomial Pj 1 o • • • o Pj t which does not vanish on p. Its degree is n — '^2 i rij i and its weights are 
given by p^ = J2i P$ • Since necessarily Pj i (p) ^ for alH = 1, . . . , I, the corresponding point 



^i" »f")=£*(4(4:> 



is contained in the convex hull displayed in the statement of the corollary. 

If P is an arbitrary covariant of degree n and weights jl^ which does not vanish on p, we can write it as a linear 
combination of monomials. If the linear combination is chosen minimally then every monomial has the same degree 
and weights as P. Since at least one of the monomials must not vanish on p, the claim follows from what we have 
proved above. □ 
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By virtue of |Corollary 2| computation of entanglement polytopes is a finite problem that can be completely 
algorithmized. Indeed, by using the relation |3"3"1 §4.2] (cf. [2]) 



C[Hf = (C[H] ® C[G/J7]) G S C[W x G/UY 



the problem of computing a set of generating covariants Pi, ... , P m is transformed into a problem of computing 
invariants for a complex reductive group, for which there exist algorithms in computational invariant theory [34 . Once 



a set of generators has been found, Corollary 2 can be used to compute the entanglement polytope both for specific 
states as well as for families of states. We use this method to compute all examples in |§ 5[ 

Finite generation also implies other desirable properties. It is clear that there are only finitely many entanglement 
polytopes, since by |Corollary 2| any entanglement polytopes is the convex hull of some subset of the finite set of points 

V = {(tlf\...,(lf ) )/n j :j = l,...,m} 
Moreover, as similarly observed in |T], the set of quantum states for which all generators are non-zero, 

m 

f]{peP(H):P 3 (p)^Q}, 

i=i 

is a finite intersection of Zariski-open sets, hence itself Zariski-open. It follows that the entanglement polytope of a 
generic quantum state is always maximal, i.e., equal to the convex hull of the finite set of all the points in V . 
Let us briefly digress to discuss this generic entanglement polytope, which we shall denote by A. Clearly, 



A = conv<{ — (jZj p 



n 



= / (/xW, . . . , juW) : 3 a covariant of degree n and weights pP^ j> 

= $(P(H)) n D x 

= {(A«...,AW) : \^=\( P ^),peP(H)}, 

the set of possible local eigenvalues of an arbitrary pure quantum state (not restricted to any particular entanglement 
class). The problem of computing this polytope is known as the one-body quantum marginal problem in quantum 
information theory and as the one-body iV-representability problem in quantum chemistry [35]. Its convexity has been 
noted in [36 - 38] , and it has been solved by combining the invariant-theoretic characterization of A with some Schubert 
calculus and geometric invariant theory |37H39) ; see also 00} for a different approach relying solely on symplectic 
geometry. In the case of N qubits, a complete solution has been obtained in |41j . 



d. Experimental Noise 

A quantum state prepared in the laboratory will always be a mixed state p and it is a priori unclear what statements 
can be inferred about its entanglement from its local eigenvalues. Here, we give two slightly different ways for leveraging 
the results discussed for pure states to the realistic noisy scenario. 

For this, we will assume that a lower bound p on the purity tr p 2 is available. One natural way of obtaining such an 
estimate is the well-known swap test ]42 which directly estimates tr p 2 using two-particle measurements on two copies 
of p. We sketch an alternative procedure which may be simpler to implement in some platforms. It is rigorous up 
to an assumption on the prevailing noise mechanism: namely that it does not increase purity. This does hold, e.g., 
for dephasing and depolarizing noise — two models applicable to the majority of experiments. Now suppose that p 
has been prepared by acting on an initial product state with a quantum operation A approximating an entangling 
unitary gate U (e.g., a spin squeezing operation). Let A' be a channel approximating U , and p' = A'(p). Under said 
assumption, we have that tr// 2 < trp 2 , i.e., the purity has decreased under the noisy "disentangling operation" A'. 
But p' is still approximately a product, so that the lower bound [43] 

N 

trp 2 >trp /2 >Y / \\* k) (p)\\l-(N-l)=P 
fc=i 
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for the global purity in terms of the local eigenvalue spectra is likely not to be too lose (it is tight for product states). 

The first way of dealing with noise is to realize that there is a pure state with fidelity (ijj\p\ij;) > p whose vector 
of local eigenvalues differs from the the vector of local eigenvalues of p by at most 2iV- v /l — p in ^-norm. To see this, 
we observe that 

p < tr pp < ||p||i HpIIoo = WpW^ = max (ip\p\ip) , 

\m\=i 

having used the matrix Holder inequality and the convention that || — || p , when applied to matrices, denotes the 
(Schatten) p-norm. Therefore, if \tp) is a pure state that maximizes the right-hand side fidelity then 

N 

53l* (k) G>) - XW(|^»||i < N\\p- MMHi < 2Ny/l - = 2N^p. (A.9) 

k=l 

The first inequality is a version of Weyl's perturbation theorem for the 1-norm 0] (11.46)], and the second inequality 
is a well-known relation between fidelity and 1-norm distance |44|, Theorem 1] . 

An alternative second approach consists in realizing that there exists a function S(p) such that if the local eigenvalues 
are more than 5{p) away from an entanglement polytope Ac, then p cannot be written as a convex combination of 
pure states |<fo) G C. More precisely, define the distance between an experimentally obtained collection of eigenvalues 
A(p) and an entanglement polytope Ac by 

N 

Then the function S(p) is characterized by the property that 

d(X(p),A c )>S(p) ptamv{\4>)(<l>\ : 0eC}. 

A simple estimate for S(p) can be derived along the lines of the previous paragraph. Indeed, assume that d (X(p), Ac) > 
S(p). Let be as above, and let \<f>) € C be arbitrary. Then 

HMM - \4M\h > - pIIi - \\p-\i>)Wh > ^ - 2v/T^ =: e(p), 



where the first step is the triangle inequality, and the second step mimics (A.9). Borrowing another estimate from |44j . 
this implies that |(-0|</))| 2 < 1 — \e{p) 2 . Now assume for the sake of reaching a contradiction that p = Pi\4>i) (4>i\ can 
be written as a convex combination of pure states \<pi) from C. Then our previous arguments imply that 



<##>=X>ilM&>| 2 <l-^tp) 2 



which contradicts the definition of whenever p > 1 — \e{p) 2 . A short calculations proves that this is certainly the 
case whenever 

8{p) > ANy/l-p. 



4. Linear Entropy of Entanglement and its Distillation 

a. Entanglement Polytopes and the Origin 

Pure quantum states p = \ip)(ip\, whose one-body reduced density matrices are maximally mixed, i.e., p^ cx 1, 
play a special role in entanglement theory. First, any such state maximizes the total uncertainty ^2j{ARj) 2 of any 
orthonormal basis of local observables Rj [T3]. Second, any entanglement monotone M(\ip)(ip\) := \P(\ip))\ defined 
in terms of a G-invariant homogeneous polynomial P attains on p its maximal value over all the states in the same 
entanglement class [T31 |19j . To see this, observe that 

M( P ) = \p(\m = \p( 9 m\ = \\ 9 i^ir \p (^^) i ^ \ p (ijfj^) 1 = M{g • p) 
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where n is the degree of P; the inequality follows from a result by Kempf and Ness which states that \\g \ip)\\ > |||V>)|| 
for all g e G if the p^ are maximally mixed |45| . 

The point in the entanglement polytopes corresponding to quantum states with maximally mixed one-body reduced 
density matrices will be called the origin and denoted by O; it satisfies oc (1, . . . , 1) for all k. Clearly O cannot be 
expressed as a proper convex combination of any two other points in the entanglement polytope, for we have arranged 



each vector of eigenvalues in weakly decreasing order. Therefore, Theorem 1 implies that the origin is contained in the 
entanglement polytope if and only if there exists a covariant with weights fp- k ' oc (1, . . . , 1). Such a covariant is of 
course nothing but a G-invariant homogeneous polynomial; hence, 

O € Ac 3 a G-invariant homogeneous polynomial which is non-zero at p. 

In particular, any entanglement monotone defined via polynomial invariants necessarily vanishes on those quantum 
states whose entanglement polytopes do not include the origin (these states are also called unstable in geometric 
invariant theory |46j : their complement being the semi-stable states). This observation has lead to the suggestion that 
unstable states should be considered "unentangled" [T^] or "not genuinely entangled" [20 , even though they might be 
entangled according to the standard notion of entanglement that we have adopted in this work. 



b. Linear Entropy of Entanglement 

The geometric picture provided by entanglement polytopes suggests another way of quantifying entanglement, which 
assigns non-zero values to all entangled states, namely the Euclidean distance of the collection of local eigenvalues. 
This quantity is precisely equal to the ^-norm distance by which the reduced density matrices differ from being 
maximally mixed. Moreover, it is by an affine transformation related to the multi-particle version of the (linear) 
entropy of entanglement [47H49] given by 

E(p) = E(\U, . . . XW) = 1-1 £||A«||f, =1-1 J2\\p^ 



k=l fe=l 



Indeed, 



I! (A«, . . . , AW) - 0||1 = 5><« - * HI = f^\\pW\\l -±=(n-J2})- NE(p). (A.10) 

k=l k k=l k \ k=l k ) 

The linear entropy of entanglement admits an operational interpretation [50J; in the case of multiple qubits, it reduces 
to the Meyer-Wallach measure of entanglement, and to the concurrence in the case of two qubits |51H53| . 

The following result then follows immediately from convexity. It naturally generalizes the properties satisfied by 
states with maximally mixed one-body reduced density matrix. 

Proposition 3. Any entanglement polytope Ac contains a unique point of minimal Euclidean distance to the origin 
O. The corresponding quantum states p maximize the linear entropy of entanglement E{p) over all states in C. 

The maximal entropy of entanglement that can be obtained from states in the closure of the entanglement class C, 
E{A C ) = max{E(p) : p e C} = max{^(A (1) , . . . A (JV) ) : (A (1) , . . . A^) e A c } 



is in view of (A. 10 1 a simple function of the Euclidean distance of the entanglement polytope to the origin and can 



thus be computed easily. 



c. Entanglement Distillation 

Given the linear entropy of entanglement E(p) as a means of quantifying multi-particle entanglement, it is natural 
to ask for a corresponding distillation procedure, i.e., a protocol for transforming the quantum state p by SLOCC 
operations to a state with maximal linear entropy of entanglement. This might only be possible asymptotically, as the 
maximum might be contained in the boundary of C. 

The function E{p) is smooth and can therefore be maximized locally by following its gradient flow in C = G ■ p. By 

B, 

dE\ p = -ld ($, . . . , pC')) \ p = -^uj FS (iX p , -) = -1 tv(X p -), (A.ll) 
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where we denote by X p the tangent vector generated at p by the infinitesimal action of the collection of reduced 
density matrices. The last identity holds since P('H) is a Kahler manifold. Therefore, the gradient of E is at every 
point p = \ip)(tp\ equal to 



. , 2 - 2d 

gradEYp = X„ = 

B yFJ N p N dt 



(e 



ip l 



mm 



t=0 



c,tp(™)^ 



(A. 12) 



which factors over the action (A. 3) of the Lie algebra of G. Therefore, any solution p t to the gradient flow equation 
remains in the entanglement class of the initial value po = p at all times t > 0. 

In mathematical terms, the Euclidean distance to the origin and hence E(p) are closely related to the norm square 
of the moment map (cf . § 3b I, and the latter is a minimally degenerate Morse function in the sense of Kirwan [3] . The 
corresponding analog of (] A .12) had already been noticed in |3J. Moreover, it was established that while the solution p t 



of the gradient flow equation might not necessarily converge to a unique limit point, E is sufficiently well-behaved so 
that E(p t ) always converges to the global maximum as t — > oo. We summarize: 

Theorem 4. By following the gradient flow of E, which at any point p is given by the infinitesimal action of the 
reduced density matrices (A. 12 1, the global maximum of the linear entropy of entanglement in the closure of the 
entanglement class of p is reached (possibly asymptotically). 



As described in the Methods Summary, |Theorem 4| is in practice implemented by approximating the gradient flow 
with finite time steps. That is, after preparing the quantum state p one measures its one-body reduced density matrices, 
re-prepares and performs local POVM measurements with Kraus operators oc e tp< ' , = y/l — (S^Y 



for sufficiently small but finite t > (cf. (A. 2)). If the outcomes of these measurements are . . . , then 

entanglement has been distilled. By successively repeating this procedure and concatenating the SLOCC operations, 
one asymptotically arrives at a quantum state with maximal entropy of entanglement. Notably, this method of 
entanglement distillation only requires local tomography and works on a single copy of the state at a time. 

From a theoretical perspective, the limit point lim t _ i . 00 p t of the gradient flow can also be seen as a normal form of 
the state p in its entanglement class (in case this limit exists). This is the point of view taken in |13| . where a similar 
numerical algorithm has been given for the special case where O is contained in the entanglement polytope. 



5. Examples 

By associating to every entanglement class its entanglement polytopes, we have obtained a finite yet systematic 
classification of multi-particle entanglement. In this section, we will illustrate their computation and application by a 
series of examples for systems of several qubits. 

Mathematically, the covariants of a multi-qubit system are in one-to-one correspondence with the covariants of 
binary multilinear forms, whose study is a prominent topic in classical invariant theory [5,4, ">"> . Before we proceed, 
we introduce some notational simplifications. It will be convenient to represent the local eigenvalues of a system of 
N qubits by the tuple (AmL, ■ ■ • , A^ x ) G [1/2, 1]^ of maximal local eigenvalues. This is of course without loss of 
information, since the sum of the two eigenvalues of any qubit density matrix, which is a 2 x 2 positive semidcfinitc 
Hermitian matrix of trace one, is equal to unity. Similarly, we may label the weights of a covariant by the tuple 
(/JmL, ■ ■ ■ , Mmax) G { \n/2~\ , . . . , n} N , since the sum of the components of each weight fl,^ has been fixed to be equal 
to the degree n of the covariant. 



a. Three Qubits 

In the case of three qubits, with Hilbert space T-L = C 2 ® C 2 (g> C 2 , there are six distinct entanglement classes [B]: 
The classes of the GHZ state |5B] and of the W state [HJ, 

ighz> = ^(ittt) + iwm i wo = ^f(ltn) + mt) + 

three classes that correspond to EPR states shared between any two of the three subsystems, 

im> = -L(itti)-mt>), \bi) = -^(im> - i;tt», |53) = -L(itit>-i4tt», 
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Covariant 


Degree 


Weight 


|GHZ) \W) 


\Bl) 


\B2) 


|B3> 


|SEP) 


/ 


1 


(1,1,1) 


X 


X 


X 


X 


X 


X 




2 


(2,1,1) 


X 


X 


X 











Hy 


2 


(1,2,1) 


X 


X 





X 








H z 


2 


(1,1,2) 


X 


X 








X 





T 


3 


(2,2,2) 


X 


X 














A 


4 


(2,2,2) 


X 


















Table I. Three-qubit covariants — labeled as in [551 p. 33] — with their degree and weight, and their vanishing behavior on the quan- 
tum states representing the six entanglement classes (x denotes non- vanishing) . The invariant A is Cayley's hyperdeterminant, 
which is closely related to the 3-tangle [141 160] . 




U 



(a) (b) (c) {d} 

Figure 5. Entanglement polytopes for three qubits: (a) GHZ polytope (entire polytope, i.e., upper and lower pyramid), (b) W 
polytope (upper pyramid), (c) three polytopes corresponding to EPR pairs shared between any two of the three parties (three 
solid edges in the interior), (d) polytope of the unentangled states (interior vertex). 



and the separable class represented by |SEP) = |ttt)- 

We shall now compute the corresponding entanglement polytopes by following the general method of covariants 
described in |§ 3 c| Using techniques crafted towards the special situation of three qubits, the same polytopes have 
already been computed in |16l 1571 155] ; the corresponding quantum marginal problem, which as we have explained 
amounts to computing the maximal entanglement polytope, has been solved in |41j . A minimal set of generators of the 
covariants of a three-qubit system (in fact, of the equivalent question for binary three- linear forms) , has been determined 



in late 19th century invariant theory [59J : There are six generators, and we have summarized their properties in Table I 
By |Corollary 2| computing the entanglement polytopes is now a mechanical task: for any quantum state representing 
the entanglement class, we merely need to collect those covariants which do not vanish on the state, and take the 
convex hull of their normalized weights. 

The resulting entanglement polytopes are illustrated in [Figure 5| They arc in one-to-one correspondence to the six 
entanglement classes described above; that is, in this particular case there is no coarse-graining. As explained before, 
one polytope is contained in the other if quantum states in the former class can be approximated arbitrarily well by 
states in the latter class. In this case, this is also a necessary condition, since there is no coarse-graining. Since the 
GHZ-class polytope is maximal, it follows that all states can be approximated arbitrarily well by states of GHZ type. 
In mathematical terms, the GHZ class is dense; this is of course well-known [BJ. Similarly, the polytope of the W class 
(upper pyramid) contains all entanglement polytopes except the GHZ one, so that by states in the W class one can 
approximate all states except those of GHZ class. 

We now illustrate our method of entanglement witnessing: 



If the point (Amax, Amax, Amix) corresponding to the collection of local eigenvalues is contained in the lower part 



of the GHZ entanglement polytope (Figure 5| (a)) 



\(!) + a (2) +A (3) <2 

'max 1 'max 1 'max ^ 

then it is by (1) not contained in any other entanglement polytope, and therefore the quantum state at hand 
must be entangled of GHZ type. 

More generally, if the point is not contained in any of the polytopes corresponding to an EPR state shared 
between two of the three particles ( Figure 5| (c), which includes (d)), i.e., if 



A£L<1 (V*= 1,2,3), 
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Figure 6. Illustration of entanglement distillation for the quantum state p = IV'XV'I discussed in § 5 a by following the gradient 
flow (curved line), the origin O = (0.5,0.5,0.5) is reached asymptotically. The corresponding quantum states have maximal 
linear entropy of entanglement. 



then by (1) the quantum state at hand must be entangled of either GHZ or W classes. These classes of states 
are the ones that possess genuine three-qubit entanglement. 

As a final example, we consider the quantum state p — where 

1 



'7 



(lttt> + im> + im) + 2|ut» 



It is easy to verify by the method of covariants that the entanglement polytope of p is full-dimensional, and by the 
above classification it follows that p is of GHZ type. However, its collection of local eigenvalues, (A,nL, A,nL, A,nL) ~ 
(0.76,0.79,0.88), is contained in the interior of the upper pyramid. As we have just discussed, the entanglement 
criterion (1) in this case only allows us to conclude that p is of either GHZ or W type. By using the entanglement 
distillation procedure described in |§ 4c| we can however transform p by into another state whose local eigenvalues are 
arbitrarily closed to the origin O = (0.5,0.5,0.5) (Figure 6 1. In this way, we arrive at quantum states which are both 



more entangled and for which our entanglement criterion (1) is maximally informative 



b. Four Qubits 



The case of three qubits was rather special, since there the entanglement polytopes represent faithfully the hierarchy 
of closures of the corresponding entanglement classes — no coarse graining takes place. In contrast, for four qubits the 
situation is the generic one: There are infinitely many entanglement classes. According to the classification of [10 , 
they can be partitioned into nine families with up to four complex continuous parameters each. Neither the families 
themselves, nor the complex parameters within these families are directly experimentally accessible. 

Here, we sketch a proof showing that the polytopal view strikes an attractive balance between reducing the complexity 
sufficiently to allow for simple experimental criteria on the one hand, and preserving a non-trivial structure on the 
other hand. Indeed, following the tradition established in |6] and continued in [10 , we may state that through the 
polytopal lense, four qubits can be entangled in seven different ways. 

We have determined all entanglement polytopes of four qubits using the general method 



§ 3 c \ applied to a minimal 



generating set of 170 covariants found in [12] . More precisely, for every family in |10) . we have analytically computed 
the covariants using a computer algebra system. Deciding whether a normalized weight ^-fi^ is included in Ac 
then amounts to solving the explicit polynomial equation Pj(p) ^ 0. The algebra system can readily perform such 
calculations analytically. 

The polytope A formed by the marginal eigenvalues of all possible pure states is the convex hull of 12 vertices. 
These can be easily described as follows: One vertex, (A 
vertices 



(i) 

max j ■ 



Amax) = (1, 1, 1, 1), corresponds to product states; six 



(0.5,0.5,1,1), (0.5,1,0.5,1), (0.5,1,1,0.5), (1,0.5,0.5,1), (1,0.5,1,0.5), (1,1,0.5,0.5) 

belong to two-partite entangled states; four vertices 

(0.5,0.5,0.5,1), (0.5,0.5,1,0.5), (0.5,1,0.5,0.5), (1,0.5,0.5,0.5) 

to three-partite entangled states, and one vertex (0.5,0.5,0.5,0.5) which is the image of a four-partite entangled state 
(not necessarily genuinely four-partite entangled — see §5c). This latter vertex is the origin as defined in §4a 
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Representative 



| GHZ) ® |t) 

\w) ® It) 

|EPR) ® |EPR) 
iEPR) ® |t) ® |t) 
It) ® It) ® It) ® It) 



Number of embeddings 



Table II. Bi-separable classes of four-qubit states. Each class corresponds to an embedding of a three-qubit entanglement 
polytope from |§ 5a| into the four-qubit polytope A in the obvious way. The right column gives the number of permutations of 
qubits that lead to distinct embeddings. 



No 


Family [TO] 


Parameters 


Vertices 


Facets 


Perms 


E(A C ) 


Vertices compared to full polytope A 


1 






12 


13 


4 


0.482143 


(0.5,0.5,0.5,0.5) replaced by (0.75,0.5,0.5,0.5) 


2 






10 


13 


4 


0.458333 


(0.5, 0.5, 0.5, 0.5), (1, 0.5, 0.5, 0.5) missing 


3 


L ai 


a = 


9 


14 


6 


0.45 


(0.5, 0.5, 0.5, 0.5), (0.5, 1, 0.5, 0.5), (0.5, 0.5, 0.5, 1) missing 


4 


Lab 3 


6 = -a =£0 


8 


16 


1 


0.5 


all three-partite entangled vertices missing 


5 


iab 3 


b = a = 


7 


9 


1 


0.375 


all three- and four-partite entangled vertices missing 


6 


L a2 b 2 


b = -a=£0 


10 


14 


6 


0.5 


(0.5, 1, 0.5, 0.5), (0.5, 0.5, 0.5, 1) missing 


7 


Gabcd 


generic 


12 


12 


1 


0.5 


n/a — this is the full polytope 



Table III. Entanglement polytopes for genuinely four-partite entangled states. The second column specifies one choice of family 
and parameters (according to the classification of [10]) which gives rise to the polytope shown (there might be further choices, 
partly because the parametrization of |10| is not always unique). "Perms" is the number of distinct polytopes one obtains when 
permuting the qubits. As defined in the main text, E(Ac) is the maximal linear entropy of entanglement in the class. The 
right-most column is a recipe for obtaining the given polytope out of the full one A. 



As in the three-qubit case, there are several lower-dimensional subpolytopes corresponding to bi-separable states. 
These are obtained by embedding the entanglement polytopes found in the previous section into the full four-qubit 
polytope in the obvious way. All possibilities are listed in |Tablc Il| 

We turn to the genuinely four-body entangled states. There are seven such polytopes, all full-dimensional. Their 



definitions and some of their properties are listed in Table III and Figure 7 The four-qubit VF-state (or Dicke state) 
|W 4 ) = ^j(ltttt) + IfN-t) + Itt-tt) + l-l-ttt}) corresponds to polytope 5. As is the case for three qubits, the polytope 
is an "upper pyramid", i.e., it is the intersection of the full polytope A with the half-space defined by 



A (l) , A (2) , A (3) , a(4) > n 



(A.13) 



Again in analogy to the three-qubit case, we can take any violation of ( |A.13 1 as an indication of "high entanglement". 
One way to make this precise is to read off Table III that violations imply that the state |?/>) can be converted into one 



with entropy of entanglement at least 0.45 (which might be much higher than E(\ip)) as obtained from the measured 
data!). The entanglement classes of the four-qubit GHZ state and of the cluster states |61_ are associated with the full 
polytope (number 7). 

There is a numerical coincidence between our findings and the ones in |20| . where also seven non-biseparable 
entanglement classes have been identified on four qubits. The two classifications are, however, not identical. Indeed, 
|20| is based purely on invariants, as opposed to the more general cewarmni-theoretic description of our polytopes. As 
mentioned in |§ 4| invariants cannot differentiate unstable entangled states from product states. Hence, the 7 classes of 
|2U] must all be semi-stable. However the same is true only for our polytopes 4, 6 and 7. In this sense, the polytope 
methods yields a finer classification for unstable vectors, whereas [2U] provides a better resolution of the stable case. 

In summary, up to permutations, there are 13 entanglement polytopes for four qubits, 7 of which belong to genuinely 
four-partite entangled classes. The numbers increase to 41 and 22, respectively, if distinct permutations are counted 
separately. 



c. Genuine k-Partite Entanglement 



In the examples so far, bi-separable states mapped onto polytopes of lower dimension. This is no longer true for 
N > 6 (for N = 6, the entanglement polytope associated with \GHZ 3 ) <g) \GHZ 3 ) is clearly 6-dimensional) . However, 
it remains true that spectral information alone can be used to witness genuine A^-qubit entanglement for any N. 



1G 















4 1, showing that our 



Figure 7. Three-dimensional cuts through the non-trivial full-dimensional entanglement polytopes for four qubits. The columns 
correspond to the first six polytopes in |Tablc III| The first coordinate is fixed in every row — to values 0.5, 0.75, and 1, respectively, 
with the remaining three variables shown. 

From the first row, one can see that only two of these entanglement polytopes (No. 4 and 6) include the maximally mixed vertex 
(0.5, 0.5, 0.5, 0.5). These reach the maximal value for the entropy of entanglement. The four polytopes miss ing the maximally 
mixed vertex cannot be distinguished from product states using polynomial entanglement monotones (c.f. 
approach is stronger in that regard. 
The final row exhibits the behavior of the class when the first particle is projected onto a generic pure state. We then recover 
a three-qubit entanglement polytope on the remaining coordinates. Polytopes 1, 3, and 6 give the full three-qubit polytope 
implying that a GHZ-type vector is generated. Polytopes 2, 4, and 5 collapse to the upper pyramid of Figure 1. Hence these 
entanglement classes produce a W-type state when the first particle is projected onto a pure state. It follows that the mixed 
3-tangle (and any other convex-roof extension of a polynomial monotone) vanishes on the mixed state generated by tracing 
out the first particle of states in these classes. This observation allows us to graphically recover some properties calculated 
algebraically in |10| : For example the vanishing 3-tangle for the class L a i 3 , a = b — (which corresponds to polytope 5). 



A general theory of genuine entanglement detection from spectral information will be presented elsewhere. Here, 
we merely give one example valid for any number N of qubits. The result is stated in the language of [62 1 (see also 
[71 151 [6UM66| and references thereirQ. Call a vector producible using k-partite entanglement if it is of the form 
W) — l"0i) ® • • * ® IVW) where every \ifji) is contained in the tensor product of at most k sites — otherwise, \ifj) is said 
to contain genuine (k + l)-partite entangled. In particular, the states that contain genuine ./V-partite entanglement are 
precisely those which are not biseparable, i.e., those which do not factorize with respect to any non-trivial (bi)partition 
of the subsystems. 

Recall first that for a system of N qubits the solution of the quantum marginal problem (i.e., the maximal 
entanglement polytope) is given by the inequalities 

A|l<£ A mL (A. 14) 

for the smallest local eigenvalues A^- n = 1 — A^L € [0, 0.5] [H|. Therefore, the polytope for the class of states that 
factorize with respect to a given partition Ax U . . . U A m — {1, . . . , N} is defined by the constraints 

*£L< E A min VieA k ,k = l,...,m (A.15) 



The reader should not be confused by the fact that some authors f ^ where semi _ stability (in the sense of pi is additionally 

use the term "genuine jV-partite entanglement" in a different sense. requ ired for "genuine" entanglement. 
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The following lemma shows that information on the local eigenvalues can serve as a witness for genuine fc-qubit 
entanglement (cf. Figure 3). In other words, certain correlations between the one-particle reduced density matrices can 
only be explained by the presence of genuine fc-partite entanglement. 

Lemma 5. Let N > 2 and k G { \N/2] , . . . , N}. For every 7 G [0, 0.5^-^] the local eigenvalues 



fA (1) 

V^min' ' 



' A min I 



1 



7 



N-r N-i' 



N - 1 



can originate only from genuinely k-partite entangled (pure) quantum states. Conversely, if k ^ N — 1 then there exists 
a realization in terms of a state producible using k-partite entanglement. 



Proof. Consider any bipartition A± U A2 = {1, . . . , N}. Without loss of generality, assume that 1 G A±. Then ( A. 15 1 
for i = k = 1 reads 



k- 1 



N ■ 



1 - 2^ 7 iv- 1 



\A, 



1 



N - 1 



\Ai\>k. 



(A.16) 



This immediately proves that any with the advertised local spectra is genuinely fc-partite entangled. For if such a 
factorizes with respect to some partition, then by (A.16 1 at least one factor must comprise at least k sites. 



To show that the local eigenvalues can indeed be realized using fc-partite entanglement, consider the bipartition 
A\ = {l,...,fc}, A2 — {k + 1, ■ ■ ■ , N}. Then (A.16 1 is satisfied. The same is true holds, trivially, for all other 



constraints for k = 1 in (A. 15 1 . The constraints for k = 2 are satisfied if and only if \A2\ ^1 <=> k ^ N — 1, which is 
true by assumption. □ 



d. N Bosonic Qubits 



The theory of entanglement polytopes is readily adapted to the case of indistinguishable particles. Here, the Hilbert 
space is no longer a tensor product, but the symmetric or antisymmetric subspace in (C d )® N , where d is the local 
dimension. Since the particles are indistinguishable, all one-body reduced density matrices are equal. Therefore, the 
entanglement polytope of an entanglement class C is in the case of indistinguishable particles defined to be 



A C = {X« : A« = A(pM), P GC} 



The reduced density matrices can be diagonalized by the coherent action of the group K = SU(<i), and the corresponding 
group of SLOCC operations is G = SL(d) [57] • Using these definitions, the theory can be developed in precisely the 
same way as above. 

We shall now illustrate this for the case of a system of N indistinguishable two-level systems with bosonic statistics, 
with Hilbert space % — Sym A '(C 2 ) C (C 2 ) (8Ar . As in the preceding examples, we may represent the vector of local 
eigenvalues of the one-body reduced density matrix by its maximum. Thus we may think of an entanglement polytope 
Ac as a line segment or interval in [0.5, 1]. 

The entanglement polytopes can be computed as above using covariants, which in this case correspond formally 
to the covariants of binary forms in mathematics; these are again well-studied and explicitly known for small N 
|54l I68| . As a convenient shortcut, we will however use a geometric argument which immediately gives the entanglement 
polytopes for arbitrary N: 

Lemma 6. For a system of N bosonic qubits, the entanglement polytopes are given by intervals [7, 1] with 

7e{i}u{^*:* = 0,l ) ...,LJV/2j 
Proof. Since the unentangled states are contained in the closure of every entanglement class, every entanglement 



polytope Ac contains the point 1. The other end point of the interval, which we denote by 7 > 0.5, is by (A. 10 1 
directly related to the maximal entropy of entanglement E(Ac), and we need to determine the possible values of 7. 
First, we observe that since the generalized GHZ state -^(If)®^ + l!) 18 ^) is locally maximally mixed, 7 — 0.5 can 
always be achieved for arbitrary N. Let us now suppose that 7 > 0.5, i.e., the spectrum of the one-body reduced 
density matrix is non-degenerate and hence contained in the interior of the set D^. Let p = |V>)(^| be a quantum state 



whose one-body reduced density matrix is a diagonal matrix with first diagonal entry 7. By (A.5), 7 can only be a 
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vertex of the entanglement polytope if the range of the differential of $ at p is orthogonal to the space of diagonal 



observables spanned by a z (otherwise, 7 could be changed infinitesimally, hence would not be a vertex). By (A.6l, 
this is the case if and only if ia zp vanishes — in other words, p is fixed by the infinitesimal action of ia z (and hence of 
any diagonal local operator) and thus \ip) is an eigenvector of the ^-operator . Such states are known as occupation 
number basis states, or Dicke states [69J. ft is easy to see that the one-body reduced density matrix of a Dicke state 
with N-f spins pointing upward and N± = N — spins pointing downward is equal to diag(iVf /N, N^/N); therefore it 
is mapped into if > . Finally, we observe that since the one-body reduced density matrix is diagonal and 
the global state fixed by diagonal local operators, every Dicke state is a maximum of the entanglement distillation 
procedure from [§"4c) It follows that 7 is indeed minimal, hence a vertex of the entanglement polytope. □ 

By associating with each of the entanglement polytopes the set of corresponding quantum states, we obtain \N/2] + 1 
families of quantum states. Each family generically consists of infinitely many entanglement classes, except for the 
unentangled one (7 = 1), and it can be characterized operationally by the maximal linear entropy of entanglement 



E(Ac) that can be achieved by states in the family. Furthermore, it follows from the proof of Lemma 6 that the Dicke 
states for 7 > 0.5 are (up to local unitaries) uniquely characterized by the property that they attain the maximal 
entropy of entanglement in their entanglement polytope. In contrast, there are generically infinitely many entanglement 
classes that contain quantum states which are locally maximally mixed (i.e., sent to the origin, corresponding to 
7 = 0.5); this follows because for N > 4 there is more than a single invariant. For example, the four qubit GHZ state 
and the Dicke state with = = 2 are both locally maximally mixed, but in a different entanglement class. 
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